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Abstract 

We suggest and demonstrate in experiment that two normal metal /insulator/ superconductor 
(NIS) tunnel junctions combined in series to form a symmetric SINIS structure can operate as an 
efficient Peltier refrigerator. Specifically, it is shown that the SINIS structure with normal-state 
junction resistances 1.0 and 1.1 kfl is capable of reaching a temperature of about 100 mK starting 
from 300 mK. We estimate the corresponding cooling power to be 1.5 pW per total junction area 
of 0.8 pm 2 at T = 300 mK. This cooling power density implies that scaling of junction area up to 
about 1 mm 2 should bring the cooling power in the /zW range. 

Introduction 

Recently it was shown [1] that in the sub-Kelvin temperature range the Peltier effect in normal 
metal /insulator/ superconductor (NIS) junctions can be used to cool electrons in the normal 
electrode of the junction below lattice temperature. The mechanism of the cooling in the NIS 
contacts is the same as that of the well-known Peltier effect in metal/semiconductor contacts - 
see, e.g., [2]. Due to the energy gap in the superconductor, electrons with higher energies (above 
the' gap) are removed from the normal metal more effectively tham those with lower energies. 

This makes the electron energy distribution sharper, thus decreasing the effective temperature 
of electron gas in the normal metal. The decrease in electron temperature demonstrated in this 
first experiment was, however, limited to about 10% of the starting temperature. There are two 
possible reasons for this. The first and most obvious, is that the resistance of the refrigerator 
junction was relatively large leading to low cooling power. The second possible reason is heat 
leakage through the SN contact used to bias the refrigerator. Nominally, an ideal SN contact 
with large electron transparency should be able to provide electric conductance without thermal 
conductance at low temperatures. However both finite temperature and finite subgap density of 
states in a superconductor lead to non-zero thermal conductance of the biasing contact which 



degrades the refrigerator performance. This poses the problem of how to bias the refrigerator 
without compromising its thermal insulation. 

The aim of the present work was to address these two problems and to show that once they 
are solved the refrigerator performance is improved dramatically. In particular, we show that 
refrigerator with two NIS junctions in series is capable of reaching temperatures of about 100 mK 
starting from 300 mK. This brings the NIS refrigerators quite close to practical applications, for 
instance, in cooling space-based infrared detectors [3]. 

The problem of large specific refrigerator junction resistance can be alleviated to some degree 
in a straightforward way by making the insulator barrier thinner. Although it is a challenging 
technological problem to push this process to its limit, we could conveniently reduce the specific 
resistance of the junctions to about 0.3 kfi x //m 2 . 

It is less obvious how to solve the second problem of heat leakage through the biasing junction. 
The solution we suggest here is to combine two NIS junctions in series to forrb a symmetric SINIS 
structure. Since the heat current in the NIS junction is a symmetric function of the bias voltage 
V, the heat flows out of the normal electrode regardless of the direction of the electric current if 
the junction is biased near the tunneling threshold, V ~ ±A/e. This means that in a symmetric 
SINIS structure we can realize the conditions when the electric current flows into the normal 
electrode through one junction and out through the other one, while the heat flows out of the 
normal electrode through both junctions. In this way the heat leakage into the normal electrode 
of the structure is minimized. The experiment with the SINIS structures described below supports 
this idea. 

Basic concepts and refrigerator structure 

We begin by briefly outlining the basic theoretical concepts concerning the heat flow in the NIS 
junctions. Under typical conditions when the transparency of the insulator barrier is small, the 
heat current P out of the normal electrode (cooling power) of an individual NIS junction is: 

(i) 

where Rj is the normal-state tunneling resistance of the barrier, fj is an equilibrium distribution 
of electrons in the jth electrode, and N(e) = 0(e 2 — A 2 ) \e\ fyjt 1 — A 2 is the density of states in 
the superconductor. From eq. (1) we can deduce several properties of the cooling power P. First 
of all, by changing the integration variable e — + — t we prove that for equal temperatures of the 
two electrodes P is indeed a symmetric function of the bias voltage, P( — V) = P(V). Plotting 



eq. (1) numerically one can see that P is maximum at the optimal bias points V ~ ±A/e. The 
optimal value of P depends on temperature and is maximum at kgT — 0.3A, when it reaches 
0.06A 2 /e 2 /?7, and decreases at lower temperatures as (kgT/ A) 3 ^ 2 [4]. Specifically, at V — A/e 
(i.e., quite close to the optimal bias voltage^ one can get from eq. (1). 



Figure 1. (a) The schematics of the SINIS refrigerator used in the measurements, and (b) an 
AFM image of the actual structure. 


Figures la and lb show, respectively, a schematic diagram of the SINIS structures studied 
in our experiments and the corresponding AFM image of the structure. Four tunnel junctions 
were fabricated around a normal metal (Cu) central electrode and four superconducting (Al) 
external electrodes. The electrodes were made with electron beam lithography using the shadow 
mask evaporation technique. The tunnel junctions were formed by oxidation in pure oxygen 
between the two metallization steps. Two junctions at the edges with larger areas were used for 
refrigeration, while the pair of smaller junctions in the middle was used as a thermometer. A 
floating measurement of voltage across the two thermometer junctions at a constant bias current 
was used to measure temperature in the same way as in the simpler one junction case [5]. 

Single-junction refrigerator 

Prior to our experiments with the SINIS refrigerator we repeated measurements in the geometry 
of Nahum et al. [1]. In our case the refrigerating junction had a resistance of Rr = 7.8 kft, the 
copper island was 10 /im long, 0.3 fi m wide and 35 nm thick. Results of the measurements of the 
electron temperature in the island as a function of refrigerator voltage V Te j T for several starting 



temperatures at V re jr = 0 are shown in Fig. 2. We see that only a few per cent refrigeration can 
be obtained, as expected. 



Figure 2. Results of the NIS single junction experiment: temperature T [mK] of the N-electrode 
versus refrigerator voltage V re f r [/iV]. Dots and solid lines show, respectively, experimental data 
and the theoretical fit including heat conductance of the biasing contact. For comparison, the 
dashed line shows the fit without the heat conductance. 

Solid lines in Fig. 2 represent a theoretical fit obtained within the standard model of electron 
energy relaxation [7-9]. Within this model, we assume that the electron-electron collision rate is 
large so that the electrons maintain an equilibrium distribution characterized by the temperature 
T, which is in general different from the lattice temperature T). The rate of energy transfer from 
electrons to phonons is then [8]: P( = EU(T S — Tf), where E is a constant which depends on 
the strength of electron-phonon coupling, and U is the volume of the island. Another element of 
the fitting process is a heat conductance k, of the biasing SN contact. (For simplicity, we neglect 
temperature dependence of k, since the temperature range of interest in Fig. 2 is not very large.) 
The value of the superconductor gap A is almost fixed by the position of the temperature dips in 
the refrigeration curves (Fig. 2) and is taken to be 155 /ieV for this sample. Solving numerically 
the equation P = P t + k, (T — T/), where P is the cooling power (1) we can calculate T as a function 
of V re f r . The fit in Fig. 2 is obtained in this way with E = 0.9 nW/K 5 /im 3 and k = 8 pW/K. For 
comparison, the dashed line shows the fit obtained for the lowest-temperature curve without /c; in 
this case E = 1.4 nW/I< 5 ^m 3 . Although there is no drastic disagreement with the data even in 




this case, we see that k improves the fit considerably. 

To see whether the value of the heat conductance k deduced from the fit in Fig. 2 is reasonable, 
we calculated the heat conductance of the SN contact with perfect electron transparency at low 
temperatures, using the method developed in [9,4]. In this approach, the heat current J in the 
contact can be written as follows: 


J = -jh- / de(e - eV)[Me - eV) - A(e)f , (c + eV) - (1 - A(e))/ 2 (e)] , (3) 

e^/t/vT J 

where Rn is the normal-state contact resistance, functions fj(e) are introduced in eq. (1), and 
A(t) is the probability of Andreev reflection from the ideal NS interface: 


| (M-(e 2 -A 2 ) 1/2 )7A 2 , M> A, 

[1, | e |< A . 


At low temperatures kgT <C A and vanishing bias voltage V eq. (3) gives for the linear heat 
conductance k: 


‘■Ww (4) 

Making use of the gap value and temperature corresponding to Fig. 2 we get from this equation 
that k = 8 pW/K corresponds to R N on the order of 100 Ohm, which is the same order of 
magnitude as in the experiment. 


Double-junction refrigerator 

Figure 3 shows our main results with the SINIS refrigerator of Fig. 1. The two refrigerating 
junctions had resistances R? = 1-0 and 1.1 kfl, respectively, and the island was 5 /im long, 0.3 fivn 
wide, and 35 nm thick. In Fig. 3 we see several refrigeration curves starting at various ambient 
temperatures at V rc j T = 0. Note that maximum cooling power is now obtained at V re j r ~ 2 A/e 
because of two junctions involved. From the position of the temperature dips we get A = 180 /ieV. 
We see that the drop in temperature is immensely improved over that of the single junction 
configuration. Solid lines in Fig. 3 show the theoretical fit which was obtained within the same 
model as for the single-junction configuration with the two modifications. We do not have heat 
conductance k this time, and we need to solve the balance equations simultaneously for the energy 
and electric current in order to determine the electric potential of the island. The best fit shown 
in Fig. 3 corresponds to E = 4 nW/K 5 /xm 3 . The fit can be classified as reasonable, although 
there are some obvious discrepancies between the data and the theory. Possible origins of these 
discrepancies include an oversimplified model of electron-phonon heat transfer on the theory side, 



and poor calibration of the thermometer toward higher temperatures on the experimental side. 
The inset in Fig. 3 shows the maximum cooling power as a function of temperature deduced 
from this fit, together with the analytical dependence obtained by summing eq. (2) over the two 
junctions. We see that the simple analytical expression (2) gives a very accurate description of 
the lower-temperature cooling power. 



Figure 3. SINIS refrigerator performance at various starting temperatures. Notations are the 
same as in Fig. 2. Dots are the experimental data, while the solid lines show the theoretical fit 
with one fitting parameter £ for all curves. The inset shows the cooling power P [pW] from the 
fits (dots), together with the analytical result from eq. (2). 

We must stress that Fig. 3 shows electron, not lattice temperature. By cooling only electron 
gas we avoid the problem of parasitic phonon heat conductance through refrigerating junctions. 
(Electron heat conduction due to tunneling plays a minor role and was accounted for in our model.) 
If we would try to cooldown the normal island as a whole, the phonon heat conduction through 
the insulator barrier would become important limiting factor also in our refrigerator. In order to 
estimate the magnitude P p u of the parasitic phonon heat flow through the tunnel junctions we 
note that neither the thin insulator layer of the junction nor the difference between parameters of 
the junction electrodes (aluminum and copper in our experiment) give rise to substantial reflec- 
tion coefficient of the large-wavelength phonons relevant at low temperatures. Therefore we can 



estimate P pk as the flow of energy associated with ballistic propagation of acoustic phonons: 

p 12071 3 c\ ’ 

where cy and c ± are longitudinal and traversal sound velocities. For aluminum at T=1 K this 
equation gives the power density P p / X ~ 0.6 nW//xm 2 . Comparing this with the cooling power of 
our refrigerator we see that the refrigerator can cool the normal island as a whole only for starting 
temperatures not greater than 0.3 K. Decrease in tunnel resistance of the refrigerating junctions 
could shift this temperature limit to few Kelvins. 


Conclusion 

In conclusion, we have shown that the nominally-symmetric SINIS structure can be used as an 
efficient Peltier refrigerator. One of the advantages of the symmetric structure is that it is easier 
to fabricate than the asymmetric single-junction configuration. Besides this, SINIS structure 
provides more efficient thermal insulation of the central electrode, which allowed us to demonstrate 
a temperature drop of about 200 mK starting from 300 mK. The achieved cooling power density 
was approximately 2 pW//im 2 , with the total power being 1.5 pW at T — 300 mK. 

The next step in the development of a practical NIS refrigerator could be further optimization 
of the refrigerator performance with respect to the resistance Rj of the insulator barrier. As we 
saw above, the cooling power of the refrigerator increases with decreasing Ry. For a fixed junction 
area this trend would continue only up to an optimal resistance, at which point the transport 
starts to be dominated by Andreev reflection [4]. The theoretical limit for the maximum cooling 
power density is on the order of 10 3 H7/^jm 2 for aluminum junctions and should be reached in the 
junctions with unrealistically low specific resistances on the order of 10~ 2 H x jx m 2 . In practice 
the limiting factors will be the technological ability to fabricate uniform tunnel barriers with high 
transparency. 
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